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Abstract 

Formulas, analogous to the Triezenberg-Zwanzig expression for the surface tension 
of a planar interface, are presented for the Tolman length, the bending rigidity, and 
the rigidity constant associated with Gaussian curvature. These expressions feature 
the Ornstein-Zernike direct correlation function and are derived from considering the 
deformation of a liquid drop in the presence of an external field. This approach is in 
line with the original analysis by Yvon in 1948. It is shown that our expressions reduce 
to previous results from density functional theory when a mean-field approximation is 
made for the direct correlation function. We stress the importance of the form of the 
external field used and show how the values of the rigidity constants depend on it. 



I. INTRODUCTION 

The most important quantity to define the properties of a liquid surface is its free 
energy or surface tension. The determination of the surface tension is therefore one 
of the key objectives in the Statistical Thermodynamical treatment of liquid surfaces 
Over the years, different routes have led to different expressions for the surface 
tension each having their own merits and limitations. 

At the start of any historical overview are mean-field expressions for the surface 
tension. At the beginning lies the squared-gradient expression derived by van der 
Waals in 1893 Q: 

oo 

a = 2m Jdz p' Q {zf , (1) 

— oo 

where po(z) is the density profile of a planar interface. The squared-gradient expression 
is a particular result derived from the more general density functional theory (DFT) 



3|-|6|, at the heart of which lies a division of the interaction potential in a hard sphere 
reference system perturbed by attractive forces, U(r) = Uh s (r) + U at t(r) . For the surface 
tension, DFT gives ]| 



c = jdzi jdf 12 U at t{r)r 2 (l - s 2 ) p' (zi) p' (z 2 ) 



(2) 



where the integration over the interparticle distance vector, ri2=r*2 — ri, is: 

2-k 1 oo 

J dr 12 = Jd(f Jds Jdr r 2 , (3) 

0-10 

which defines s = cos(#i2) and z 2 = z\ + rs. A gradient expansion of the density in 
the expression for the surface tension in Eq. (j2j) yields the van der Waals expression in 
Eq.(Tj]) with the identification, m— — (1/12) Jdfi2 r 2 £/" a tt (?") • 

A second approach results in the so-called virial expression for the surface tension, 
first derived by Kirkwood and Buff in 1948 ?]]: 



(4) 



(2) 

where p (zi, z 2 , r) is the pair density correlation function of a planar interface. Even 
though this expression is derived by Kirkwood and Buff using a mechanical approach, 
by way of the evaluation of the components of the pressure tensor in the interfacial 
region it is an exact equation, the only assumption being made is that of pair- 
wise additivity of the molecular interactions. The theoretical evaluation of the surface 
tension using the Kirkwood-Buff expression requires knowledge of the pair correlation 
function in the interfacial region, which is hard to access theoretically but can be 
determined to great accuracy in computer simulations [8j, |9|. 

A third approach, which also provides an expression for the surface tension that is 
exact, results in the so-called Triezenberg-Zwanzig (TZ) expression 
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111: 



a 



dz x \dfi2 r 2 (l - s 2 ) p' (zi) p' (z 2 ) C (zi, z 2 , r) 



(5) 



where Cq{z\, z 2 , r) is now the direct correlation function 12] of the planar interface. 

In this article, we are interested in deriving expressions for the surface tension of a 
curved surface. In particular, we investigate the expansion of the surface tension o~(R) 
to second order in the (inverse) radius of curvature R as described by the curvature 



coefficients 5 (the Tolman length) 13] , k (the rigidity constant associated with bending 
or bending rigidity) and k (the rigidity constant associated with Gaussian curvature). 
For a spherical and cylindrical surface, this expansion has the form [fjj]: 

, m 25a (2k + k) . , . (n . 

a s (R) = a - — + R2 +... (sphere) (6) 

/ t-,\ 5a k . , 

a c (R) = a -^ + ^ + --- (cylinder) 

The virial expressions and mean-field expressions (DFT) for 5, k and k have been 
derived some 20 years ago 15 but until now, the corresponding expressions in 
terms of the direct correlation function are lacking. Our goal in this article is to fill 
this void and to derive Triezenberg-Zwanzig-like expressions for the Tolman length 
5, and the rigidity constants k and k. To achieve this, we follow in the footsteps of 
the pioneering work by Jim Henderson and coworkers who addressed precisely this 



problem in three papers in the early 80's 18M2CM] . They derived an expression for the 



Tolman length in terms of the direct correlation function, which is consistent with the 
expression presented here. Furthermore, they pointed to difficulties in extending the 
analysis to second order in the expansion in the curvature, i.e. to the order of the 
rigidity constants k and k, which find their origin in the fact that the thermodynamic 
conditions used to impose a certain curvature of the interface then play a role. 

The article is organized as follows. In Section [TT] we address the Thermodynam- 
ics involved in deriving expressions for the rigidity constants in general. Two distinct 
thermodynamic routes are formulated which correspond to the two original derivations 
of the Triezenberg-Zwanzig expression for the surface tension of the planar interface: 
(i), the fluctuation route, used by Triezenberg and Zwanzig themselves in 1972, in 
which one considers the second order change in free energy due to surface fluctua- 
tions; (ii), the equilibrium route, used by Lovett et al. [211 ]. following an approach 
suggested by Yvon in 1948 [10], in which one investigates the deformation of a planar 
interface induced by an external field. These two approaches, which result in the same 
expression for the surface tension of the planar interface, are discussed in Section II III 
and extended to curved surfaces. The full calculation of the Tolman length 5 and the 
rigidity constants k and k is presented in Section IIVI with details of the calculations 
left to the Appendices. We end with a discussion of results. 
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II. FLUCTUATION ROUTE VERSUS EQUILIBRIUM ROUTE 



Historically, the Triezenberg-Zwanzig expression for the surface tension has been de- 
rived in two distinct ways: (i), by the analysis of the second order change in free 
energy due to surface fluctuations [ll|; (ii), by applying an external field to deform 



an initia 



ously [22 



l y p lanar interface 10|, l2l| . These two routes, which we have termed previ- 



23| as the fluctuation route and the equilibrium route, are connected 



to the (thermodynamic) path considered to deform the surface and to induce a certain 
interfacial curvature. 

In the fluctuation route, the curvature of the interface naturally comes about by 
the presence of thermal fluctuations. The bulk region is unaffected and the chemical 
potential remains that at two-phase coexistence, fi = fi coex - This route is in many 
senses equivalent to a route in which one considers the presence of a non-uniform 
external field that is non-zero in the interfacial region only. The "fluctuation" is then 
envisioned as being induced by the non-uniform external field. 

In the equilibrium route, the interface is curved by changing the value of the chem- 
ical potential to a value off- coexistence at fixed temperature (fi > yU CO ex)- One then 
considers the surface tension c r (R) of a spherically (or cylindrically) shaped liquid 
droplet with radius R in metastable equilibrium with a bulk vapour. In effect, a uni- 
form external field is added of the form V ext {r) = — A/i, where Afi = fi — /i coex - Since 
the uniform external field acts throughout the system, it also affects the bulk densities 
far from the interfacial region. 

or the planar surface tension, 
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-[20|] that this may not be 



These two distinct routes lead to the same expression 
but it was already noted by Henderson and coworkers 
the case for the radius dependent surface tension c r (R), and in particular for the 
coefficients in an expansion in 1/R of the surface tension beyond the Tolman length. 
Their concerns were connected to the fact that different external fields may lead to 
the same interfacial curvature and that it is not a priori clear that the expression 
for the radius dependent surface tension derived using one external field equals that 
of another. Although the two routes do lead to the same expression for the surface 
tension, it turns out that this indeed is not the case for the rigidity constants. This was 
demonstrated previously in the context of DFT expressions for the bending rigidity 



k 



23 



241 ] : the fluctuation route was investigated by supposing the presence of an 
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external field chosen proportional to the derivative of the density profile to ensure 
that it acts in the interfacial region only 23|. It was shown that the value obtained 



for the bending rigidity then differs in magnitude and scaling behaviour from that 



obtained using the equilibrium route 



24J (see Appendix [A} . 



That the thermodynamic conditions imposed to induce a certain interfacial curva- 
ture have bearing on the rigidity constants, may also be deduced from the fact that 
the rigidity constants depend on the density profile p%(z), which describes the leading 
order change in the distribution of molecules when the interface is curved. The fact 
that the density profiles p\{z) are different in the fluctuation and equilibrium route, 
follows directly from the observation that in the equilibrium route, p\{z) is not equal 
to zero in the bulk regions since the chemical potential is also different (/i 7^ /i CO ex), 
while it is zero in the fluctuation route since fluctuations of the interface do not affect 
the thermodynamic state of the bulk regions. 

With this important observation in mind, we consider, in the next section the usual 
derivations of the Triezenberg-Zwanzig expression for the surface tension of a planar 
interface and discuss how they are extended to determine the Tolman length 5 and 
the rigidity constants k and k. 



III. TWO DERIVATIONS OF THE TRIEZENBERG-ZWANZIG EXPRES- 
SION FOR THE SURFACE TENSION 

To understand how the direct correlation function enters in the Triezenberg-Zwanzig 
expression for the surface tension, we need to address the Statistical Thermodynamical 
definition of the direct correlation function 121 ]. 



In the fluctuation route considered first, the direct correlation function describes 
the (second order) response of the free energy due to density fluctuations. This route 
was used in the original derivation by Triezenberg and Zwanzig [Uj in 1972. In the 
equilibrium route, the direct correlation function describes how a small change in the 
local density can be interpreted as a change in the local external field (second Yvon 
equation). This approach, inspired by the work of Yvon in 1948 [loj], was used by 



Lovett et al. [2l| to rederive the TZ expression for the surface tension of the planar 
interface by considering the external field necessary to bend an initially flat surface into 
a spherical droplet. This procedure was then extended by Henderson and Schofield 



5 



19[ to spherical droplets in order to have access to the full radius dependent surface 
tension. 



A. Fluctuation route 



In the fluctuation route, one investigates the change in free energy due to density 
fluctuations around a planar interface. The direct correlation function describes the 
second order change in the free energy in response to such a density fluctuation: 

k B T 



An 



dn df 12 C (z 1} z 2 , r) 8p{fx) 5p(f 2 ) . 



(7) 



The density fluctuation can be expressed in terms of a height function h(r\\) that 



describes the shift in the location of the surface 



23 



25|: 



5p(r) = p(r) - p {z) = -p' {z) h{f\\) - 



v 2 Mni) + ... , 



where it is assumed that | V/i|<^C 1 and where pi(z) is identified as the correction to 
the density profile due to the curvature of the interface. In general, pi{z) describes 
the rearrangement of molecules that results when the interface is curved, but its form 
still depends on the way the curvature is brought upon. 

The surface tension and bending rigidity are thermodynamically defined as the 
change in free energy due to height fluctuations [26] : 



AQ = - 
2 



dq 



aq 2 + kq A + ...}h(q)h(-q). 



(9) 



Inserting Eq.(JH]) into Eq.([7]) and comparing the result to Eq.Q, we retrieve the TZ 
expression for the surface tension in Eq.(j5]). A further result is that for the bending 
rigidity one finds 25]: 

k B T " 



k 



4 



k B T 



4 
k B T 



64 



dzi Jdfi2 pi(zi)pi(z 2 ) Coizt, z 2 ,r) 

) 

OO 

Jdzi Jdf l2 r 2 (l - s 2 ) p' (z 1 )p 1 (z 2 ) C (zi, z 2 ,r) 
dzi df 12 r 4 (l - s 2 ) 2 p' (zi)p' (z 2 ) C (zi, z 2 ,r) 



(10) 



This formal expression for k was first derived by Parry and Boulter 25J . Even though 
this expression is exact, for the explicit evaluation of it, one still requires some way of 
determining the density profile p\{z). Furthermore, expressions for 5 and k cannot be 
derived using this route. 



B. Equilibrium route 



To understand the way the direct correlation function enters the expression for the 
surface tension in the equilibrium route, we need to discuss the two Yvon equations 
which are in the heart of any Statistical Thermodynamics treatment of surfaces l[ Q, 
0]. In the first Yvon equation, the total correlation function, G(rxJ$i)i relates a small 
change in the local external field to a change in the local density % |5[ : 

1 



5p(n) 



df 2 G(n,f 2 ) SV ext (f 2 ) 



(11) 



where G(fi,r 2 ) is related to the pair density as: 



G(fi,f 2 ) = p {2) (n,r 2 ) -p(fi)p(f 2 ) +p(fi)S(f 2 -fi) 



(12) 



The inverse of the total correlation function defines the direct correlation function 12 j 
and it describes in the second Yvon equation how a small^change in the local density 
can be attributed to a change in the local external field 



^Kxt(^i) = -k B T Jdf 2 C{f 1 ,f 2 ) Sp(f 2 ) . 



(13) 



By considering translational deformations, the second Yvon equation appears in the 
following form: 

ViV^Fi) = -k B T Jdf 2 C(n,f 2 ) V 2 p(f 2 ) , (14) 
which, in planar geometry, reduces to: 



KxtOi) = ~ k sT Jdr 2 Cq{z u z 2 ,t) p' (z 2 ) . 



(15) 



As a next step, we consider the influence of the external field on the Laplace pressure 
difference Ap across a spherical or cylindrical droplet with radius R. The generalized 
Laplace equation for the pressure difference across a surface with (constant) curvatures 
J=l/R 1 + l/R 2 and K = 1/(R 1 R 2 ) is given by |27|: 

k 



Ap = aJ-25aK--J 3 + 2kJK + ... 

2 



(16) 



For a spherically or cylindrically shaped surface, this expansion takes the form: 

(sphere) (17) 



2a 25a 



Ap 



a 



k 



R 2i? 3 



(cylinder) 



18) 



where the dots represent terms of order 0(1/R 4 ), which indicates that the term pro- 
portional to 1/i? 3 is absent in the expansion of the spherical interface. It is important 



to mention that these expressions are derived with the equimolar radius |28j chosen 
as the radius of the spherically (or cylindrically) shaped surface, i.e. R = R e , where 
R e for a spherical liquid droplet is defined by: 

4vr Jdr r 2 [p,(r) "/>«] = y Rl (pi ~ Pv) ■ (19) 


Next, we consider the consequences of mechanical balance in the presence of an 
external field. The condition of mechanical balance can be expressed in terms of the 
pressure tensor [l|: 

V-pV) = -p(f)W cx t(r). (20) 
In spherical symmetry, the expression for mechanical balance leads to: 

Pnir) = ~[p T (r) -p N (r) ] - p s (r) V^{r) , (21) 
r 

which integrated from inside the liquid to the vapour phase gives 

Ap = Jdr -[ p N (r) - px(r) ] + Jdr p s {r) V^ t {r) . (22) 


following the derivation by Yvon [10], this expression was applied by Lovett et al. 

21j to rederive the TZ expression for the surface tension by considering the situation 
in which an initially flat surface, in the absence of an external field, is bend into 
a spherical droplet by an infinitesimally small external field SV ex t(r) in such a way 
that the bulk regions are unaffected. More generally, one may consider an already 
curved spherical droplet, whose curvature is increased or decreased (R — > R + SR) 
by the infinitesimally small external field [ijij]. In either case, the change in pressure 
difference is zero, 5(Ap) = 0, and the two terms on the right-hand-side of Eq. (1221) must 
balance. Since we have that in the absence of the external field 

OO 

r, 2 r , , , , , 2cr 28a , . 

Ap = Jdr -[ p N (r) - p T (r) ] = — - — + ..., (23) 

this implies that 

R R 2 







*(^-^ + --)= ldr P ' s (r)5V ext (r), (2 J) 
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where we have integrated the last term in Eq. (|22p by parts. The second Yvon equation 
in Eq.f lT3|) is then used to rewrite the right-hand-side as: 



2a 45a 



R 2 R 3 



\ CXJ 

+ ..A5R=k B T jdn jdr 2 C s {r u r 2 ,r) P ' a {ri)5p 8 {r 2 ). (25) 



The leading order term in the expansion in 1/R of the right-hand-side in this equation 



leads to the TZ expression for a 



2l|, w : 



an expression for the Tolman length 



rile the next to leading order term leads to 



18 



19) . Here, we are interested in the second 



order terms beyond the leading term. Even though the coefficient of the 1/R 4 term in 
Eq.()25l) happens to be zero, its absence shall provide information for the determination 
of the combination of rigidity constants 2k + k. 

In order to obtain an expression for the bending rigidity k separately, the analysis 
is also carried out for a cylindrically shaped interface. In cylindrical geometry, the 
consideration of mechanical balance than ultimately leads to 



a 3k 



R 2 2R 4 



LXJ 

SR-wfafacJn.wMMi^. (26) 



For the cylindrical interface, the 1/R 4 term is related to the bending rigidity k and 
it is the 1/R 3 term that vanishes. In the next section, Eqs. fl25]) and (1261) are used to 
carry out a systematic expansion in 1/R to derive explicit expressions for the Tolman 
length 5, and the rigidity constants k and k. 



IV. EQUILIBRIUM ROUTE TO DETERMINE 6, k, AND k 

We consider the situation outlined in the previous section: a spherical liquid droplet 
with radius R is deformed by an infinitesimally small external field into another spher- 
ical droplet with radius R^-R + SR without affecting the bulk regions. To determine 
the resulting change in density profile, we use the fact that we can write the density 
as p s (r) = p s (r — R;R). This means that when R is changed, the density profile is 
affected in two ways: (i), the variable z = r — R is shifted and (ii), the density profile 
itself depends on R. This means that 

5 Ps (r) = p s (r - R-5R;R + 5R) - p s (r - R; R) 

= -p' s (r) 5R + p 3 {r;R + SR) - p s (r; R) . (27) 
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When we expand p s {r) in 1/R, as in Eq. flA3p . this gives 
5p s {r) = 



SR. 



{21 



R 2 R 3 

This expression is inserted into Eq. fl25j) and expanded to second order beyond the 
leading order term. Details of the analysis for spherical droplets are outlined in Ap- 
pendix [Bj details of the analysis for cylindrical droplets in Appendix O For the Tolman 
length, the analysis of the spherical interface in Appendix [B] leads to Eq. (1B9j) : 



5a 



Jdz 1 Jdr 12 r 2 (l - s 2 ) [p' s (z x )p' s {z 2 ) C s (z u z 2 , r) ] x 



(29) 



where the subscript 1 refers to the leading order correction in an expansion in 1/R (see 
Eq. (IA3j) ). The evaluation of the Tolman length using this expression therefore requires 
knowledge of the way the density profile and direct correlation function depend on 
the radius of curvature. 

The expression for the Tolman length in Eg . 0291) can also be derived from the 

nri 

expression for a s (R) presented by Henderson and coworkers in refs. |18l. |19|]: 
k B T 



CT S (R) 



dri / df 12 r 2 - (n - r 2 ) 2 ] p' s {r 1 )p' s (r 2 ) C s (r 1 ,r 2 , r) 



(30) 



which is constructed to be correct to order 0(1/R) [3]. An explicit expansion in 1/R 
of the term in square brackets gives: 



r 2 - (r x - r 2 ) 



r 2 (l 



R ! 



(31) 



Inserting this expansion in Eq.( l30l and comparing with the expansion for a s (R) in 
Eq.Q results in the expression for the Tolman length in Eq. fl29l) . It should be men- 
tioned that starting with the Henderson expression in Eq. (l30l . an expression for the 
Tolman length very similar to Eq.f l29|) was derived by Barrett in 1999 29|. The only 
difference is that in the Barrett expression for the Tolman length, C s> x(zi, z 2 , r) in 
Eq.( )29|) is rewritten in terms of the triplet direct correlation function of the planar 
interface 29 1. 

When the expression for the Tolman length in Eq.( l29|) is combined with the result 
of the analysis of the cylindrical interface in Appendix [Cl we find as an alternative 
expression for the Tolman length: 

k B T 



5a 



(32) 
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It is reminded that in this expression the location of the z = plane corresponds to the 
equimolar surface. The advantage of this expression for the Tolman length is that it 
can be evaluated from the properties of the planar interface only. In Appendix |A] we 
verify that both expressions for the Tolman length in Eqs. (T29]) and f|32|) are consistent 
with known DFT expressions 13, I2J \ ■ 

The vanishing of the 1/R 4 term in Eq. (|25p in the analysis for the spherical interface, 
leads to Eq. (IBlOp in Appendix [B] This gives the following two expressions for the 
combination 2k + k: 

k T °° 

2k + k = -- jL- fdz x Jdf 12 r 2 (l - s 2 ) [p' s (zi)p' s (z 2 ) C s (zi, z 2 , r) ] 2 (33) 

—00 

2k B T 



3 

—00 



knT 



knT 



dzi Jdf 12 rs p' (zi)p Sj2 (z 2 ) C (zi, z 2 ,r) 
dz x df 12 rsp 1 (z 2 ) [p' s {zi)C s (z 1 ,z 2: r)) 1 



-00 
00 

Idzt Jdf 12 r 2 (l - s 2 ) z\ p' Q {zi)p' Q {z 2 ) C (z 1 ,z 2 ,r) 
—00 

k T °° 

+~^~ j dzi J df ^ 12 r4 ( X ~ Po( z i)p'o( z 2) C (z u z 2 ,r) , 



and 



00 

2k + k = — |- ldz l jdf 12 r 2 {l-s 2 )z 1 [p'^zMz^C^z^r)], (31) 



—00 
00 



k B T 

3 

—00 

00 



dzi df 12 rsz! p' (zi)pi(z 2 ) C (zi,z 2 ,r) 



k B T 
' 12 



dz! df 12 r (1 - 3s ) p' G (zi) pi(z 2 ) C (zi, z 2 , r) 



-00 
00 



+ ~W J dZl J d ^ 12 ~~ s4 "> Po( z i)Po( z 2) C o(zi, z 2,r) . 

—00 

The analysis for the 1/R 4 term in Eq.( )26|) for the cylindrical interface, leads to 
Eq. (IC9j) in Appendix O This gives the following two expressions for the bending 
rigidity k: 

k T °° 

k = — |— jdzi df 12 r 2 (l - s 2 ) sinV [p' c {z\)p' c {z 2 ) C c (z h z 2 ,(p,r)] 2 (35) 
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and 



4A: B T 
3 

k B T 
6 

k B T 
4 

k B T 
' 64 



dzi Jdf 12 rs p' (zi)p Ci2 (z 2 ) C (zi, z 2 ,r) 
rfzi df 12 rsp!(z 2 ) [p' s (z 1 )C s (z 1 ,z 2 ,r)] 1 



— oo 
oo 



-oo 
oo 

oo 
oo 

dzi df 12 r 4 (l - s 2 )(l + 3s 2 )po(z 1 )po(2 2 )Co(^i,2 2 ,^) 



k B T 
12 



oo 

2/1 „2 



6 



dz 1 Jdf 12 r (1-s )^i [p / s (zi)p / s (z 2 )C s (zi,z 2 ,0]i (36) 

dZl /df 12 rs^i Po(2i)pi(22) Co (21, 3s, r) 



— oo 
oo 



k T 

+ /^i jdfu r 2 (l -3s 2 )p' (2;i)pi(2; 2 )Co(2;i,2; 2 ,r) 



-oo 
oo 



k T 

1 — |— /^i /dr 12 ^ 2 (1 - s 2 )z 2 p' (2i)po(z 2 )Co(^i,2; 2 ,r) 



-oo 
oo 



y^i2 r4 (! - •s 2 ) 2 p' (2;i)p[ ) (2; 2 )Co(2i,22,r) . 

— oo 

Combining Eqs. (l34p and (1361) then leads to the following expression for k: 
k T °° 

k = — ^— jdzi df 12 r 2 (l - s 2 ) zf p' (zi)p' (z 2 ) C (zi, z 2l r) (37) 



4 

L,/ Uzi ldf 12 r 4 (l - s 2 )(l - 5s 2 )po(2:i)po(z 2 ) C (^i, z 2 , r) . 



— oo 

fc R T 



96 

— oo 

The expression for the Tolman length in Eq. fl3"2"|) and the expressions for the rigidity 
constants in Eqs.()36l) and (|37j) are the final results of this article. 



V. DISCUSSION 

For the surface tension of the planar interface, one may distinguish three type of 
microscopic expressions: 

• Mean-field or DFT expressions of which the van der Waals squared-gradient 
expression is an example. These are approximate in nature but are widely used 
since they are relatively easy to evaluate (numerically). 
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The Kirkwood-Buff virial expression, which is exact and widely used to deter- 
mine the surface tension in computer simulations. 

The Triezenberg-Zwanzig expression in terms of the direct correlation function. 
Although the evaluation of the surface tension using this expression requires 
access to the direct correlation function, it has the advantage that it is an exact 
expression and the short-ranged nature of the direct correlation function has 
proven to be helpful to settle fundamental issues such as the influence of capillary 
waves on the structure and tension of a planar interface in the limit of vanishing 



gravitational field 3CH32|. 



There exists a great level of consistency between these three type of expressions for 
the surface tension. Although it is quite elaborate to show the equivalence of the 

3311 , it is 



Triezenberg-Zwanzig and Kirkwood-Buff expression for the surface tension 
quite more simple to show how they are both consistent with mean-field expressions. 
Ever since the introduction of the coefficients S, k and k that appear in the curvature 



expansion of the surface tension by Tolman [13| and then later by Helfrich [ijj] , it has 
been the goal to formulate the same type of microscopic expressions and to achieve the 
same level of consistency as it exists for the surface tension of the planar interface. The 
road to the realization of this goal started with the formulation of squared- gradient 
expressions, first for the Tolman length by Fisher and Wort is in 19 84 34J ] and then 



for the rigidity constants by Gompper and Zschocke [3_5f| and by us 16j in the early 



90 's. These squared-gradient expressions were then shown to be embedded in more 



general DFT expressions in 1998 



13, 



24|. By that time, also the virial expressions 



for 5, k and k, analogous to the Kirkwood-Buff formula for the surface tension, were 
formulated 15] and it was verified that they reduce to the DFT expressions when a 



mean-field approximation for the pair density is made. 

Important progress was also made in the formulation of Triezenberg-Zwanzig-like 
expressions in terms of the direct correlation function. Even before the work by Fisher 



and Wortis in 1984, Henderson and coworkers [18M20I] derived an expression for o~(R) 
correct to first order in the expansion in 1/R. Although it was not shown explicitly at 
the time, the expression for the Tolman length in Eq. fl29l can be extracted from their 
analysis, even though the reduction to the more simple expression in Eq.( j32l) could not 



yet be made. Further progress was then made by Parry and Boulter [251 ]. who extended 



13 



the original analysis by Triezenberg-Zwanzig and derived the analogous expression for 
the bending rigidity in Eq. ffTO]) . To understand the status of this expression, it is 
important to realize that the very nature of the rigidity constants k and k as the 
change in free energy due to the curvature of the interface, immediately leads to the 
following two rather subtle issues: 

• The curvature of the surface, as described by the total and Gaussian curvatures 
J and K, is not uniquely defined simply because the allocation of the position of 
the dividing surface cannot be made unambiguously. Always, some procedure to 
locate the dividing surface either by an integral constraint (equimolar surface) or 
by a crossing constraint has to be chosen. For the surface tension (and Tolman 
length), this observation bears no consequences but the rigidity constants are 
intrinsically linked to some choice for the location of the dividing surface. For 
example, the expressions for k and k in Eqs. fl36|) and fl37|) are derived using the 
equimolar surface as dividing surface (Eq. (1A12l) ). 

• The system's response to curvature, and therefore the value of the rigidity con- 
stants, depends on the way this curvature is induced. Most straightforwardly 
it is induced by uniformly changing the value of the chemical potential {equilib- 
rium route) - the mean-field and virial expressions are all derived in this way - 
but one could also envision the application of a non-uniform external field {fluc- 
tuation route). The corresponding values of the rig idit y constants then differ as 



demonstrated previously in the context of DFT [22|,[36|]. The implication is that 
the equivalence of the equilibrium route and the fluctuation route, as it exists for 
the derivation of the Triezenberg-Zwanzig expression for the surface tension of 
the planar interface, does not exist for the derivation of Triezenberg-Zwanzig-like 
expression for the rigidity constants. 



The expression for the bending rigidity by Parry and Boulter [25| is derived in the 
fluctuation route, which leaves the derivation of Triezenberg-Zwanzig-like expressions 
for the Tolman length S and the rigidity constants k and k in the equilibrium route 
as the final piece of the puzzle. In this article, the expressions sought after are finally 



derived by extending the analysis by Henderson and coworkers [18l-l20j . For the Tolman 
length, our final result is the expression in Eq.( )32|) . For the rigidity constants our final 
expressions are presented in Eqs.f l36|) and ([3j 
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The Triezenberg-Zwanzig-like expressions for the Tolman length and rigidity con- 
stants possess, like the Triezenberg-Zwanzig expression for the surface tension of the 
planar interface itself, certain advantages with respect to other approaches. Foremost, 
the TZ-like expressions are exact expressions which means that they are valid beyond 
the mean-field approximation. This is not merely a quantitative issue since it is well- 
known that the presence of capillary waves are not fully incorporated in mean-field 
theory as it fails, for instance, to capture the divergence of the interfacial width in 
the limit of vanishing gravitational field {4, 30]. A further advantage of the TZ-like 



expressions is that the direct correlation function featured is short-ranged even in the 
situation that the capillary length diverges and the range of the pair density corre- 
lation function becomes infinite. This is instrumental in addressing the influence of 
capillary waves on the surface tension of non-planar interfaces. For the surface tension 
of the planar interface it is well-established that capillary waves reduce the mean-field 
value for the surface tension by as much as 20%, but the same analysis for the Tolman 
length or the rigidity constants has not been carried out thus far. 
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Appendix A: Density functional theory expressions 

In this appendix, we present the results of density functional theory for 5, k and 



k (details of the analysis can be found in ref. 24j) and show how the Triezenberg- 
Zwanzig-like expressions in terms of the direct correlation function reduce to the DFT 
expressions when a mean-field approximation is made. 

The expression for the (grand) free energy in DFT is based on the division into a 
hard-sphere reference system plus attractive forces described by an interaction poten- 
tial ?7 a tt (t) ■ It is the following functional of the density 3_-f|: 

Q[p] = fdf [ / hs (p) - w(r) } + - fdn fdf 12 U M (r) p{n)p{r 2 ) , (Al) 



2 

where /h s (p) is the free energy density of the hard-sphere reference system with uniform 
density p. This functional is based on a local density approximation for the hard- 
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sphere reference fluid, but more sophisticated approaches exist using weighted-density 
6| or fundamental measure theory [37|. 

The Euler-Lagrange equation that minimizes the above free energy is given by: 

fM =»- Jdf 12 U att (r) p(r 2 ) . (A2) 

To derive expressions for the curvature coefficients 5, k and k, one needs to consider 
the expansion of the density and free energy in Eq.f lAlj) to second order in 1/R for 
spherically and cylindrically shaped liquid droplets. For instance, the expansion of 
the density profile of the spherical droplet reads: 

Ps{r) = po{z) + ps,i{z) + PsM + ■■■ ■ ( A3 ) 

The leading order correction to the density profile of the spherical interface is twice 
that of the cylindrical interface, so it is convenient to define pi(z) = p s i(z) = 2 p c tl (z) 
where z = r — R. 

The coefficients in the curvature expansion of the density are determined from the 
curvature expansion of the Euler-Lagrange equation in Eq. (1A2I) . The result is that 
the (planar) density profile po(z) is determined from solving: 



Pcoex = /hs(Po) + jdr 12 U^r) p {z 2 ) , (A4) 
and pi(z) follows from solving: 

2 

Pi = fM Pi(*i) + fern ^att(r) [ Pl (z 2 ) + ^(1 - s 2 ) p' (z 2 ) ] , (A5) 



where p\ = 2a/ Ap |l6, 36]. For the evaluation of the curvature coefficients it turns 
out to be sufficient to determine the density profiles po(z) and Pi(z) only. 

Expressions for the curvature coefficients are now derived from the expansion of 
the free energy in Eq.f lAll) to second order in 1/R and by comparing the result to the 
expansions in Eq.([6]) . This leads to the following two equivalent expressions for the 
Tolman length [171 ]: 



1 X 



5a = ^ jdzi jdr 12 C/ att (r)r 2 (l - s 2 ) p'^z^p^fa) , (A6) 



and 



-oo 



oo oo 



Sa = - Jdzi Jdf 12 f/ att (r)r 2 (l - s 2 ) z 1 p' Q (zi) p' G {z 2 ) - y \dz z p' (z) , (A7) 



-oo — oo 
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The second order term in the expansion of the free energy for the spherical interface, 
eads to the following expression for the combination of the rigidity constants, 2k + k 



2k + k 



dz x df 12 U att (r)r (1 - s ) p' (zi)pi(z 2 ) 



-oo 
oo 



- 7 /(fej jdf 12 U att (r)r 2 (l - s 2 ) z\ // (>i) p'o(z 2 ) 

— oo 
oo 

+ - jdz l jdr 12 U. Att {r)r\l-s*)p' Q {z l )p' Q (z 2 ) 



+ dz 



yzp'i(z) +pi z 2 p' {z) + p Sj2 zp' (z) 



(A8 



where p s>2 = —aAp 1 /(Ap) 2 — 25a / Ap [16|, [36| with Api = p X £ — p\ >v . An alterna- 
tive expression that contains no reference to the chemical potential is also derived in 
ref. |2^: 



2k + k = ~ jd Zl jdr l2 C/ att (r) r 2 (l - s 2 ) p'^p'^Zi) 



— oo 
oo 



-- jdzx \dr l2 C/ att (r) r 2 (l - s 2 ) p'^z,) p[(z 2 ) 



-oo 
oo 



-- Jdzx Jdf 12 U att (r) r 2 (l - s 2 ) z\ p' (zi)p' (z 2 ) 

— oo 

1 00 

+ - Jdzt Jdf 12 U att (r) r 4 (l - s 4 ) p' Q {z x ) p' Q {z 2 ) . 



(A9) 



To derive an expression for the bending rigidity k, an expansion of the free energy 
to second order in 1/R is made for the cylindrical interface. For the Tolman length the 
expressions in Eqs. (1A6I) and (1A7I) are recovered and one finds for the bending rigidity 



k = - jdzildr 12 U altt {r)r (1 - s ) Pq{zi)p 1 {z 2 ) 



(A10) 



+ ^ jd Zl ldf 12 U att (r)r\l- s 2 ) 2 p' ( Zl )p' (z 2 ) 



oo 

+ [dz 



^zp' l (z) + ^z 2 p' (z) + 2p c>2 zp' (z) 
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where p C)2 = —a Ap x / (2 Ap)' 
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As an alternative expression for k that contains 



no reference to the chemical potential is 



oo 

k = - J dzi J dr 12 U &t t(r)r 2 {I - s 2 ) p'q(zi)p' c 2 (z 2 ) 



{All) 



•i j dr.. jdf l2 I (;• ) r 2 { 1 - s 2 ) ,/(; )//, ( : 2 ) 



-oo 



-~ jdz x jdf 12 U att (r) r 2 (l - s 2 ) 2 z\ pf Q (z x ) p' (z 2 ) 



-oo 
oo 



+ __ j dzi j d f l2 Uatt ( r ) r 4 (1 _ s 2 )(1 + 3s2) p i^ Zl)p i^ Z2 ) . 



These expressions are derived in ref. [24| without reference to any particular choice 
for the location of the dividing surface. In the present analysis, it is important that 
the radius R is defined according to the equimolar radius in Eq. ffTUj) . Expansion 
of Eq.( fT9l) to first order in l/R then leads to the following two conditions for the 
expanded profiles po(z) and pi{z) \l\ : 

oo oo oo 

Jdz z p'q(z) = , and Jdz z p[(z) = — Jdz z 2 p' (z) . (A12) 

— oo — oo — oo 

By making explicit use of this property for po(z), Eq. ( 1A7I) may be rewritten as: 



^=7 / dz i i d ?i2 U att (r) r 2 (l - s 2 ) z x p' (zx) p' {z 2 ) . 



(A13) 



Furthermore, using the properties for po(z) and pi(z) in Eq. (lA12j ). the expression for 
2k + k in Eq.( 1A8l) is rewritten as: 



2k + k = ~ \dz x jdr l2 U att (r) r 2 (l - s 2 ) {z x + z 2 ) p' { Zl ) p' x {z 2 ) (A14) 



while Eq. (IA10p is rewritten as: 



k = -- d Zl df 12 U att (r) r 2 (l - s 2 ) {z x + z 2 ) p' (z 1 )p' 1 (z 2 ) (A15) 



-oo 

oo 



+ - dzi df 12 U att {r) r (1 - s ) z x p' (z x ) p' (z 2 ) 



+ __ j dzi f d f 12 f/ att ( r )r 4 (l - s 2 ) 2 p'^p'^z,) . 
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Next, we verify that the expressions for 5, k and 2k + k as presented in Section 
IIVI reduce to the DFT expressions listed above. To show the equivalence of the 
Triezenberg-Zwanzig and DFT expression for the surface tension of the planar inter- 
face, one approximates the direct correlation function by the attractive part of the 
interaction potential 0, l^]: 



k B TC (z 1 ,z 2 ,r) 



(A16) 



Inserting this approximation into the TZ expression for the surface tension in Eq.(j5]) 
then immediately leads to the DFT expression in Eq.(|2]). 

Eq. flA16l) is easily generalized to non-planar geometries by replacing Cq and po by 
the corresponding non-planar direct correlation function and density profile. Since all 
the curvature expressions feature a factor r 2 times the direct correlation, we may dis- 
regard the delta-function contribution in Eq. flA16j) . Using the fact that the interaction 
potential is independent of R, we have the following replacements: 



k B T W a C a ] x 
^bT [p' s p' s C s } 2 
kBT [p' c p' c C c } 2 



p'^U^r), (A17) 
[f/ (ziWz*)+p{(zi)f/ (z2)] U att (r), 
Pq(ziK, 2 (*2) + p^iK^a) + ^(siKteO ] U att (r) , 

Po(^l)Pc, 2 (^2) + pi(^l)pi(>2)/4 + ^,2(^1)^0(^2) I EWO • 



Another property that is helpful in verifying that the expressions in Section HVl reduce 
to the DFT expressions, is integration by parts with respect to the parameter s. It is 
easily verified that: 



x 1 2 

Jdsrsp x (z 2 ) = Jds r —{l-s 2 )p'^{z 2 ), 

-1 -1 
1 1 

Jds(l-3s 2 )p x (z 2 ) = - Jdsrs(l-s 2 )p' x (z 2 ) 



(A18) 



With the help of Eqs. (lA17jl and ( 1A18I) . we may now show that: (i), For the Tolman 
length 5: Eq.(J22J reduces to Eq.^SBJ and Eq.flBS]) reduces to Eq.( fAT3l) : (ii), for the 
combination 2k + k: Eq. (l5Hi) reduces to Eq. f1A9[) and Eq. flM]) reduces to Eq. (lA14P ; 
(iii), for the bending rigidity k: Eq.( l35|) reduces to Eq. (1A11[ ) and Eq. (l36l) reduces to 
Eq.(TAl5D. 
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Appendix B: Equilibrium route — Spherical droplet 



In this appendix, we insert the expression for Sp s (r) in Eq. (!28|) into Eq. (l25l) and expand 
to order 0(1/ R 4 ). Three terms result that we shall investigate separately: 



2a Aoa 



R 2 R 3 



-k B T Jdr 1 Jdr 2 C s (r h r 2 ,r) p'^n) p' s (r 2 ) (Bl) 



dn df 2 C s (r h r 2 ,r) p' s (ri) pi(r 2 ) 



dri df 2 C s (n,r 2 ,r) p' s (ri) p St2 (r 2 ) . 



R 2 

o 

2k B T 



The second Yvon-equation in Eq. (|T4|) in spherical symmetry reads: 

Kxt( r i) n = fdr 2 C s (r u r 2 , r) p' s {r 2 ) h . (B2) 

Considering the component in the f\ direction and considering the limit of a vanishing 
external field, we have that: 

(ri + sr) 



= Jdf 2 C s (n, r 2 , r) p' s (r 2 ) 



r 2 



(B3) 



where we have used that fi ■ r 2 = (r\ + sr)/r 2 . This is conveniently rewritten as: 

(ri + sr 



dr 2 C s (ri, r 2 , r) p' s (r 2 ) = df 2 C s (n, r 2 , r) p' s (r 2 ) 



r 2 



The term in square brackets is expanded in the inverse radius to give: 



(B4) 



(ri + sr) 
r 2 



2R 2 
2 



^ ' _2 , n . r /-i - 2 , 



I - -{ Zl + sr) + ^ Ui + 2^ lS r - — (1 - 5s') j + ... 



(B5) 



where it is reminded that z\=r\ — R. For the second and third term in Eq.( 1Bll) . we 
can use that 



dr x df 2 C s (ri, r 2 , r) p' a (r x ) p x (r 2 ) 



dn df 2 C s (r u r 2 , r) p' s (r 2 ) p x (ri) 



(ri + sr) 



r 2 



r 2 



(B6) 



The first term in the square brackets results after making use of the 1 <-> 2 symmetry 
of the direct correlation function. The second term in the square brackets can be 
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added since Eg. (IB4j) indicates that it is equal to zero. This addition is convenient 
since the sum of these two terms is now of order 0(1/ R): 



r\ [r\ + sr) 



r\ 



2sr 2z\sr 
+ 



:i-7s 2 ) + 



(B7) 



r 2 R R 2 2R 2 

An important subtlety in the expansion of the direct correlation is the fact that the 
second argument in C s (t\ , r 2 , r) is expanded around r 2 =r± + rs. The way this is done 



is explained in ref. 



15| where the same mathematical technique is used for the pair 



density. The result is that: 



C 8 (ri,r 2 ,r) 



1 + 



r 2 (l 



d 



r (1 — s 



2\2 



2R rds 
r 2 (l — s 2 ) z 2 d 



+ 



+ . . . 



d 2 



8R 2 r 2 ds 2 



C s (r x ,rx + rs,r) . 



(B8) 



2R 2 r ds 

A systematic expansion in 1/R of the right-hand-side in Eq.f lBlj) can now be made. 
After some algebra (using Eq. (1D4I) in Appendix [D| and by comparing with the left- 
hand-side of Eq.( 1Bll) . we obtain the TZ expression for a in Eq.fjSJ) and find for the 
Tolman length the following expression: 



5a 



JdziJ dr 12 r 2 (l - s 2 ) [p' s (z 1 )p' s (z2)C s (z 1: z 2: r)\ l . 



(B9) 



From the l/R 4 term in the expansion of the right-hand-side of Eq.( 1Bl[) . we have: 



0= Jdz 1 Jdf 12 r 2 (l-s 2 ) [pt.{zi)pt.{z 2 )C a {z 1 ,Z2, 



2k B T 



jdz-i_jdr 12 rs pi(z 2 ) [p' s (zx)C s (zx,z 2 ,r)] l 



-co 
oo 



Jdzijdfu rs p' Q (z x )p Sj2 (z 2 ) C (z 1} z 2 ,r) 



— oo 
oo 



>2k + k (BIO) 



^ Jd Zl Jdf 12 r 2 (l - s^zlp'^p'^C^z^r) 



+ 



+ 



+ 



48 



-oo 
oo 



Jdzijdr 12 r 4 (l - s 4 ) ^0(^1)^0(^2) C (zi, z 2 , r) 



—00 
00 



Jdz!jdf 12 r(l -s 2 )zi [p' s (zi)p' s (z2)C s (z 1 ,z 2 ,r)} 1 



—00 
00 



fdzxfdru rs z x p' Q (z^)p x (z 2 ) C (z 1} z 2 , r) 



-00 
00 



Jd Zl Jdr 12 r 2 (l -3s 2 ) p' ( Zl ) Pl (z 2 )C ( Zl , z 2 ,r) 



> - (2k + k) 



-00 
00 



^ fd Zl fdr 12 r 4 (l - s 4 ) p (^i)Po(^) C ( Zl , z 2 , r) 



The brackets indicate how the vanishing of this term may expressed as the sum of 
two separate expressions for 2k + k. This division may appear rather arbitrary (the 
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fifth and ninth term cancel, for instance), but is constructed such that the separate 
contributions reduce to the DFT expressions for 2k + k as outlined in Appendix |A] 
This means that in principle a term could be added and subtracted that has the 
property that it is zero when a mean-field approximation is made, such as a term that 
contains C s ,i(zi, z 2 , r) only. Since none of the terms in Eq. (1B10l ) have this property, 
the presence of such a term is hard to imagine, however. 



Appendix C: Equilibrium route Cylindrical droplet 

In this appendix, we carry out the same analysis as in Appendix [B] but now for the 
cylindrical interface. Expanding p c (r) in 1/R, as in Eq. fl28|) . gives 



$Pc{r) 



P>(r) + P -M + M 4 



SR. 



(CI) 



2R 2 R 3 

Inserting this expression for Sp c (r) into Eq.( [26l) and expanding to order C(l/i? 4 ), we 



now have: 



a 3k 



R 2 2R A 



00 

-k B T Jdri Jdf 2 C c (ri,r 2 ,(f,r) p' c (ri) p' c (r 2 ) 



(C2) 



dn df 2 C c (n,r 2 ,(p,r) p' c (ri) pi(r 2 ) 



dn df 2 C c (n, r 2 , (p, r) p' c {rx) p c , 2 (r 2 ) 



k B T 
2R 2 



2k B T 



In the absence of an external field, the second Yvon-equation in Eq. (H4|) in cylindrical 
geometry gives that: 

(ri + sr) 



jdf 2 C c (n, r 2 , (f, r) p' c (r 2 ) = Jdf 2 C c (ri, r 2 , (p, r) p' c (r 2 ) 



1 - 



r 2 



(C3) 



where we have used that V 2 pc{ r 2) — p' c i r 2) r 2 and that ri • r 2 — (r x + sr) / r 2 . The term 
in square brackets can be expanded in the inverse radius to give: 



(ri + sr) 



r 2 



2R 2 



[1 - s 2 ) sinV 



(C4) 



x 



1 ~ R^ Zl + sr ^ + i 2 i^ 1 + 2ziSr ~ T^ 1 ~ ^ sil)2(p + r2s2 ) + 



For the second and third term in Eq. (IC2l) . we can use that 

X) 
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dr x df 2 C c (r u r 2 , <p, r) p' c (r 2 ) p x (n) 



r i iji + sr) 
r 2 r 2 



(C5) 



The first term in the square brackets results after making use of the 1 4-)- 2 symmetry 
of the direct correlation function. The second term in the square brackets can be 
added since Eq. flC3l) indicates that it is equal to zero. This addition is convenient 
since the sum of these two terms is now of order 0(1/R): 
r i ( r i + sr) 



sr Zisr r 2 s 2 
r 2 r 2 ~~R + ~R 2_ + 2R 2 + "' ' 

Analogous to Eq. flB8j) . C c (ri, r 2 , (f, r) is expanded around r 2 = r\ + rs: 

,2\ • 2,„ a ^4/i n 2\2 „• 4,„ j2 



C e (r 1 ,r 2 ,(p,r) 



r 2 (l — s 2 ) sm 2 (f d r 4 (l — s 2 ) 2 sinV d 2 
+ 2R Tds + 8R 2 r^s 2 

r 2 (l — s 2 ) z 2 sin 2 <£> d 



(C6) 



(C7) 



C c (r 1 ,r 2 + rs,(p,r) . 



2R 2 r ds 

Again, a systematic expansion in 1/R of the right-hand-side in Eq. flC2j) can now be 
made. After some algebra and by comparing to the left-hand-side of Eq. (IC2l) . we now 
obtain from the condition of the vanishing of the first order term that: 



= ML - r .i. /-./- ,.2/-i ..2 



fdz 1 fdr 12 r 2 (l - s 2 ) z\ p' (zi) p' (z 2 ) C (z 1: z 2 , r) 



.k B T 



-oo 
oo 



JdziJdfu r 2 (l - s 2 ) [p' s (zi)p' s (z2) C a (z u z 2 ,r)} 1 



5a (C8) 
-5a 



where we have used the fact that to leading order C C) i(zi, z 2 , <p, r) = | C S) i(zi, z 2 , r). 
The 1/R 4 term in Eq.( lC2l) now supplies us with an expression for the bending rigidity: 

oo 

k= *f- Jd Zl Jdr 12 r 2 (l - s 2 ) sm 2 <p [f/ e {z 1 )f/ c (z 2 )C c (z 1 ,z 2 ,(p,r)] i 



4fc B T 



Jdz 1 fdf 12 rs p'oiz^p^z^ C Q {z 1 , z 2 , r) 



+ 



+ 



k B T 



k B T 



— oo 
oo 



Jdzijdfu rspi(z 2 ) [p' s {z x ) C s {z 1 , z 2 , r) ] t 



-oo 
oo 



!dzxjdr 12 r 2 (l - s 2 ) z\ f/ (zi) p' (z 2 ) C (z 1} z 2 , r) 



-oo 
oo 



fd Zl fdr 12 r\l - s 2 )(l + 3s 2 ) p' G { Zl ) p' (z 2 ) C (z u z 2 , r) 



-oo 
oo 



k (C9) 



-Y Jdz.Jdn.^il-s^z, [p' s {z 1 )p' s {z 2 )C s {z 1 ,Z2,r)} l 



jdzijdr 12 rs z x p'^z x )p x (z 2 ) C {zi, z 2 , r) 



-oo 
oo 



+ 



+ 



M 

12 



— oo 
oo 



Jdzijdr 12 r 2 (l - 3s 2 ) p' (zi) pt(z 2 ) C (zi, z 2 , r) 



\ 2k 



-oo 
oo 



Jdz 1 Jdf 12 r 2 (l - s 2 ) z\ p' ( Zl )p' (z 2 ) C (z u z 2 , r) 



— oo 



+ ML fd Zl fdf 12 r 4 (l - OVo(*iK(*2)0>(*i,32,r) 
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Again, the brackets indicate how the bending rigidity may expressed as the sum of two 
separate expressions. This division is constructed such that the separate contributions 
both reduce to the DFT expressions in Appendix [A] but, again, a term could be added 
and subtracted that has the property that it is zero when a mean-field approximation 
is made. 



Appendix D: Properties of direct correlation function 

In this appendix, we investigate some properties of the direct correlation function 
of the planar interface, Cq(z\, z 2 , r), which turn out to be useful in rewriting the 
expression for the Tolman length in Appendix [B] 

In a uniform external field, the second Yvon equation in Eq. ffl5|) reduces to: 



= k B T J dr 12 p' (z 2 ) C (zi, z 2 , r) . (Dl) 

The second Yvon equation in Eq. ({T3]) applied to the situation in which a planar 
interface is curved by shifting the chemical potential off-coexistence leads to: 

2 

t i 1 = k B T Jdf 12 [ Pl (z 2 ) + s 2 ) p' (z 2 )}C (z 1 ,z 2 ,r) , (D2) 

Eq.( 1Dll) multiplied by z\ p\{z,\) and integrated over z\ gives: 

oo 

= k B T Jdzi J df 12 (z x + rs) (J Q {z\) px{z 2 ) C (z u z 2 , r) . (D3) 

— oo 

Eq. (ID2l) multiplied by z\ p'q{z\) and integrated over z\ gives: 

Idzx / dr X2 rs p' (zi)pi(z 2 ) C (z!, z 2 ,r) (D4) 



1 x 



2 

— oo 

where Eq. flD3l) has been used. 



dz x \ df 12 r (1 - a ) Zi p' Q {zi) p' {z 2 ) C (z 1 , z 2 , r) , 
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